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Abstract. The solar photospheric Fe abundance has been de- 
termined using realistic ab initio 3D, time-dependent, hydrody- 
namical model atmospheres. The study is based on the excel- 
lent agreement between the predicted and observed line pro- 
files directly rather than equivalent width, since the intrinsic 
Doppler broadening from the convective motions and oscil- 
lations provide the necessary non-thermal broadening. Thus, 
three of the four hotly debated parameters (equivalent widths, 
microturbulence and damping enhancement factors) in the cen- 
ter of the recent solar Fe abundance dispute regarding Fe I lines 
no longer enter the analysis, leaving the transition probabilities 
as the main uncertainty. Both Fe I (using the samples of lines of 
both the Oxford and Kiel studies) and Fe II lines have been in- 
vestigated, which give consistent results: logeFei = 7.44±0.05 
and log eFcii = 7.45 ± 0.10. Also the wings of strong Fe I lines 
return consistent abundances, logeFcii = 7.42 ± 0.03, but due 
to the uncertainties inherent in analyses of strong lines we give 
this determination lower weight than the results from weak and 
intermediate strong lines. In view of the recent slight downward 
revision of the meteoritic Fe abundance log ep c — 7.46 ± 0.01, 
the agreement between the meteoritic and photospheric values 
is very good, thus appearingly settling the debate over the pho- 
tospheric Fe abundance from Fe I lines. 

Key words: Convection - Hydrodynamics - Line: formation - 
Sun: abundances - Sun: granulation - Sun: photosphere 



1. Introduction 

The solar iron abundance is of fundamental importance as it 
provides the standard to which all other elemental abundances 
in stars are compared. Furthermore, since Fe is the dominating 
contributor to the total line-blanketing and a significant elec- 
tron donor for late-type stars such as the Sun, the exact value 
of the Fe abundance influences the overall photospheric struc- 
ture. Thus Fe indirectly affects the emergent spectrum and the 
derived abundances for other elements as well. In spite of its 



significance and after many previous investigations, the solar 
Fe content is still, astonishingly enough, debated on the level 
of 0.2 dex (Blackwell et al. 1995a,b; Holweger et al. 1995). 
Not even the basic reasons for this large discrepancy have been 
properly understood, though it is commonly blamed on differ- 
ences in adopted parameters for the analysis (g/-values, equiv- 
alent widths, collisional damping parameters, microturbulent 
velocities) as well as subtle differences in computer codes 
(Kostik et al. 1996). Such dissonance even for the Sun natu- 
rally rises concern regarding derived stellar abundances with 
claimed accuracies of 0.05 dex or less. 

Recently, with the advent of accurate ^/-values for weak 
Fe II lines (Heise & Kock 1990; Holweger et al. 1990; Biemont 
et al. 1991; Hannaford et al. 1992; Raassen & Uylings 1998; 
Schnabel et al. 1999) and improved treatment of the colli- 
sional damping of Fel lines (Milford et al. 1994; Anstee et 
al. 1997), there seems to be some convergence towards find- 
ing consistency between the photospheric and the meteoritic 
Fe abundances (Grevesse & Sauval 1998, 1999). The stud- 
ies by the Oxford group (Blackwell et al. 1995a,b and refer- 
ences therein), however, stand out with their distinguished high 
value of log efoI = 7.64 [] rather than the current best esti- 
mate of log epc = 7.50 for the meteoritic abundance (Grevesse 
& Sauval 1998, but see Asplund 2000, hereafter Paper III) as 
determined from carbonaceous chondrites of type 1 (CI chon- 
drites). 

It is sobering to remember that all of the above-mentioned 
investigations rely on several approximations and assumptions 
not necessarily justified in the case of the Sun. Traditional 
abundance analysis of stars are based on one-dimensional (ID), 
theoretical model atmospheres constructed under the assump- 
tions of plane-parallel geometry (or spherical geometry for 
stars with extended atmospheres), hydrostatic equilibrium (or 
steady state stellar winds for hot stars), flux constancy and with 
the convective energy transport computed through the mix- 
ing length theory (Bohm-Vitense 1958) or some close relative 
thereof (e.g. Canuto & Mazzitelli 1991), with all their limita- 
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tions and free parameters. For late-type stars the simplifying 
assumption of LTE is also normally adopted (cf. Gustafsson 
& J0rgensen 1994 for a review of stellar modelling of late- 
type stars). For the Sun, information of the emergent spectrum, 
e.g. details of the limb-darkening, may be utilized to construct 
a semi-empirical model atmosphere, such as the widely used 
Holweger-Muller (1974) model (the assumptions of ID, plane- 
parallel geometry, hydrostatic equilibrium and LTE have here, 
however, been retained). 

A closer inspection of the solar photosphere reveals that 
none of these assumptions are strictly correct: the solar sur- 
face is dominated by the granulation pattern reflecting the con- 
vection zone deeper inside, which results in an evolving in- 
homogeneous surface structure with prominent velocity fields 
between warm upflows (granules) and cool downflows (inter- 
granular lanes) with very different temperature gradients (e.g. 
Stein & Nordlund 1998). Of course there are also regions with 
significantly enhanced magnetic field strengths, which may in- 
fluence the emergent spectrum. It is therefore not surprising 
that none of the available ID model atmospheres, including 
the Holweger-Muller (1974) model, satisfactory predict simul- 
taneously all of the various observational diagnostics (limb- 
darkening, flux distribution, H-lines etc) even for the Sun (e.g. 
Blackwell et al. 1995a; Allende Prieto et al. 1998) Naturally, 
different species are affected differently by the granulation and 
its heterogeneous nature. Lines from the dominant ionization 
stage (Fe II in the case of the Sun) will be less influenced by 
the details of the atmospheric structure while lines from other 
ionization stages are more sensitive. Furthermore, lines from 
minority species will in general be more susceptible to depar- 
tures from LTE, which can be expected to be more pronounced 
in inhomogeneous atmospheres compared with ID model at- 
mospheres (cf. discussion in Kostik et al. 1996). 

Furthermore, abundance analyses normally proceed with 
additional assumptions when synthesizing the spectrum. In or- 
der to approximately account for the photospheric velocity 
fields and produce the needed extra line broadening, both a 
microturbulent velocity £ tU rb - supposedly representing small- 
scale velocities - as well a macroturbulent velocity - reflect- 
ing large-scale motions not present in the model atmospheres - 
are applied to the spectral synthesis. The exact shapes of these 
additional broadening recipes to be convolved with the syn- 
thetic spectrum also remain poorly understood (cf. Gray 1992). 
Finally, the treatment of collisional line broadening normally 
stems from the approach by Unsold (1955), enhanced by an 
ad-hoc factor to account for the still lacking amount of broad- 
ening for strong lines. With recent quantum mechanical cal- 
culations (Anstee & O'Mara 1991, 1995; Barklem & O'Mara 
1997; Barklem et al. 1998) the introduction of the unknown 
damping enhancement factor may no longer be necessary how- 
ever, at least not for lines of neutral species. 

There is therefore no doubt that traditional abundance 
determination are built on somewhat shaky grounds, which 
need to be verified by more detailed calculations. With the 
recent progress in ab-initio numerical multi-dimensional hy- 
drodynamical simulations of surface convection of stars (e.g. 



Stein & Nordlund 1989, 1998; Nordlund & Dravins 1990; 
Atroshchenko & Gadun 1994; Freytag et al. 1996; Kim & Chan 
1998; Asplund et al. 1999a; Ludwig et al. 1999; Trampedach 
et al. 1999) there is fortunately an alternative to classical 
model atmospheres. Such inhomogeneous atmospheres self- 
consistently calculate the convective energy transport and the 
velocity and temperature structures, making the concepts of 
mixing length parameters, microturbulent and macroturbulent 
velocities obsolete. The high-degree of realism of these con- 
vection simulations is supported by the fact that they suc- 
cessfully reproduce the solar granulation pattern and statistics 
(Stein & Nordlund 1989, 1998), helioseismological constraints 
such as p-mode frequencies and depth of the convection zone 
(Rosenthal et al. 1999), and spectroscopic diagnostics such as 
flux-distribution, limb-darkening and detailed line profiles and 
asymmetries, even on an absolute wavelength scale (Asplund 
et al. 1999b; Asplund et al. 2000b, hereafter Paper I). In the 
present paper we apply such granulation simulations of the Sun 
to the problem of the solar Fe abundance, utilizing both weak 
lines and the wings of strong lines. In order to minimize the 
impact of possible departures from LTE both lines of neutral 
and ionized Fe have been investigated. By fitting the line pro- 
files the uncertainties introduced with equivalent widths can be 
avoided. Furthermore, whenever possible the improved colli- 
sional broadening treatment of Anstee & O'Mara (1991) has 
been used. 

2. 3D model atmospheres and spectral line calculations 

The procedure for calculating the spectral line transfer is the 
same as in Paper I and therefore only a short summary will be 
given here. For additional information on the details of the con- 
vection simulations and the 3D spectral synthesis, the reader is 
referred to Paper I. 

Realistic ab-initio numerical hydrodynamical simulations 
of the solar surface convection have been performed and used 
as 3D, time-dependent, inhomogeneous model atmospheres 
with a self-consistent description of the convective flow and 
temperature structure in the photosphere. A state-of-the-art 
equation-of-state (Mihalas et al. 1988) has been used together 
with the 3D equation of radiative transfer which included the 
effects of line-blanketing (Nordlund 1982) with up-to-date con- 
tinuous (Gustafsson et al. 1975 with subsequent updates) and 
line opacities (Kurucz 1993). The original simulation has a res- 
olution of 200 x 200 x 82, which was interpolated to a grid with 
dimension 50 x 50 x 82 to ease the computational burden in the 
spectral line calculations. Simultaneously the vertical resolu- 
tion was improved by only extending down to depths of about 
700 km compared with the initial 2.9 Mm. Various test ensured 
that this procedure had no effect on the resulting profiles. The 
convection simulation used for the spectral synthesis here and 
in Paper I covered about 50min on the Sun. For the present 
purposes the time coverage is sufficient to obtain properly spa- 
tially and temporally averaged line profiles, as verified by test 
calculations; even intervals as short as lOmin result in abun- 
dances within 0.02 dex of the estimates using the whole time- 
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sequence. The resulting effective temperature is very close to 
the nominal solar value, T c s = 5767 ± 21 K, while the adopted 
surface gravity was \ogg = 4.437 [cgs]. For the computa- 
tion of background continuous opacities and equation-of-state, 
a standard solar chemical composition was used (Grevesse & 
Sauval 1998). In particular the assumed He abundance (10.93) 
was consistent with the helioseismological evidence (e.g Basu 
1998; Grevesse & Sauval 1998) though the exact value is of no 
practical importance for the present investigation. 

In the present investigation only intensity spectra at solar 
disk center (p, — 1.0) were considered, which have been calcu- 
lated for every column of the snapshots, before spatial and tem- 
poral averaging and normalization. The assumption of LTE in 
the ionization and excitation balances and for the source func- 
tion (S v = B v ) have been made throughout in the line transfer 
calculations. The line profiles were computed for 141 velocities 
around the laboratory wavelength with a interval of 0.2 (weak 
and intermediate strong lines) or 1.5-2.0kms _1 (strong lines); 
one additional point was computed without consideration of the 
line to estimate the continuum intensity necessary for the nor- 
malization. All lines were calculated with three different abun- 
dances (logeFc = 7.30, 7.50 and 7.70) from which the final 
profile with the correct line strength was interpolated from a 
X 2 -analysis of the whole profile in a similar fashion to the study 
of Nissen et al. (2000); test calculations ensured that the abun- 
dance step was sufficiently small not to introduce any signifi- 
cant errors in the derived abundances (Alog e-pc -C 0.005 dex). 



3. Line data and observed solar spectrum 

The accuracy of the final results naturally depend not only on 
the degree of realism of the model atmosphere but also on the 
quality of the necessary input atomic data. The choice of tran- 
sition probabilities for the lines is a delicate matter, with many 
recent discussions of pros and cons in the literature. The gf- 
values for Fe I lines of the Oxford group (Blackwell et al. 1995a 
and references therein), the Hannover-Kiel workers (Holweger 
et al. 1995 and references therein) and O' Brian et al. (1991) 
are of excellent internal consistency and all agree to within 
0.03 dex on average (Lambert et al. 1996). The last source has 
larger quoted uncertainties in general, which is also evident in 
the significantly larger scatter in the derived abundances. For 
completeness we have included both remaining samples in the 
analysis; the overlap of lines with quoted equivalent widths 
W\ < 10 pm is limited to nine lines, which differ by 0.028 dex 
on average. However, given the in general good reputation of 
furnace measurements and the fact that concern regarding the 
quality of the (^/-values from the Hannover- Kiel group recently 
has been voiced (Kostiketal. 1996; Anstee et al. 1997), we tend 
to give the results obtained with the Oxford transition probabil- 
ities greater weight. Also the scatter in the derived abundances 
is slightly larger when adopting the Hannover-Kiel ^/-values. 
When deriving the Fe abundance from the wings of strong Fe I 
lines, we have been guided by the quality measures quoted by 
Anstee et al. (1997) and selected the most suitable lines (in total 



Table 3. The adopted line data for the Fe II lines 



Wavelength a 




log<?/ b 


log7 rad a 


w x b 


log e Fc 


[nm] 


[eV] 






[pm] 




457.63334 


2.844 


-2.94 


8.612 


6.80 


7.42 


462.05129 


2.828 


-3.21 


8.615 


5.40 


7.35 


465.69762 


2.891 


-3.59 


8.612 


3.80 


7.40 


523.46243 


3.221 


-2.23 


8.487 


8.92 


7.49 


526.48042 


3.230 


-3.25 


8.614 


4.74 


7.63 


541.40717 


3.221 


-3.50 


8.615 


2.76 


7.38 


552.51168 


3.267 


-3.95 


8.615 


1.27 


7.35 


562.74892 


3.387 


-4.10 


8.487 


0.86 


7.49 


643.26757 


2.891 


-3.50 


8.462 


4.34 


7.38 


651.60716 


2.891 


-3.38 


8.464 


5.75 


7.52 


722.23923 


3.889 


-3.36 


8.617 


2.00 


7.60 


722.44790 


3.889 


-3.28 


8.617 


2.07 


7.55 


744.93305 


3.889 


-3.09 


8.612 


1.95 


7.28 


751.58309 


3.903 


-3.44 


8.612 


1.49 


7.49 


771.17205 


3.903 


-2.47 


8.615 


5.06 


7.41 



a From lohansson (1998, private communication) and the VALD 
data base (Kupka et al. 1999) 

b From Hannaford et al. 1992. Note that W\ is only listed here to al- 
low easy identification in Fig. |5| and is not used for the abundance 
determinations 



14 lines). The ^/-values for these lines are taken from Black- 
well et al. (1995a) and O'Brian et al. (1991). 

For the Fell lines there are five recent sources for gf- 
values: Heise & Kock(1990), Hannaford et al. (1992), Biemont 
et al. (1991), Raassen & Uylings (1998) and Schnabel et al. 
(1999) of which the first two and the last are based on experi- 
mental data while the remaining two have been obtained from 
semi-empirical calculations. Again, the variations between the 
different compilations are relatively small on average, though 
occasionally large on a line-by-line comparison. We tend to 
view the theoretical calculations with some balanced scepti- 
cism due to the noticably larger scatter in derived abundances 
when selecting the g/-values by Biemont et al.; for the nine 
lines in common between all five sources, the standard devia- 
tion increases from 0.07 dex for the values by Hannaford et al. 
and Heise & Kock to 0.13 dex when using Biemont et al.'s pre- 
dictions (cf. also discussions in Hannaford et al. 1992 and Bell 
et al. 1994). Furthermore, with the data from Biemont et al. the 
derived Fe abundances show a distinct trend with wavelength, 
suggesting a problem in the calculations. With the more recent 
calculations by Raassen & Uylings the scatter is improved to 
a comparable level to the measured ^/-values, but the abso- 
lute scale is clearly offset (— O.lOdex relative to Hannaford et 
al.) compared with the other four sources, and therefore we are 
hesitant to adopt these calculations here. An investigation of 
the reason for these differences seems worthwhile (cf. Grevesse 
& Sauval 1999). The final choice between the remaining three 
compilations is somewhat arbitrary, but we have opted for the 
measurements by Hannaford et al. (1992) as it includes two 
additional lines (15 lines in total) and the absolute scale of 
their gf- values is in between the other two sources; adopting 
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instead the values by Heise & Kock (1990) and Schnabel et 
al. (1999) would change the derived abundance by +0.04dex 
and -0.02 dex, respectively, but leave the line scatter essentially 
unaltered. We note that the more recent lifetime measurements 
of Schnabel et al. have slightly smaller claimed uncertainties 
than in Hannaford et al.. However, since the independent ex- 
periments of Guo et al. (1992) support the measured lifetimes 
of Hannaford et al., we will here retain the Hannaford et al. gf- 
values but keep in mind that the derived Fe II abundances may 
be overestimated with 0.02 dex on average. 

For the collisional broadening from hydrogen atoms the 
quantum mechanical calculations developed by Anstee & 
O'Mara (1991) have been applied for the Fe I lines for transi- 
tions between levels of type s-p, p-s, p-d, d-p, d-f, and f-d. The 
broadening cross-sections and their dependence on tempera- 
ture have been kindly provided by Barklem (1999, private com- 
munication) from line-by-line calculations (the data has subse- 
quently been incorporated into the VALD database, Barklem et 
al. 2000). For a few lines not individually computed, the nec- 
essary data was obtained from interpolation in tables provided 
by Anstee & O'Mara (1995), Barklem & O'Mara (1997) and 
Barklem et al. (1998). The contribution from collisions with 
helium atoms have been included by assuming that the cross- 
sections scale with the polarizability of the perturbing atom 
in the same way as in the van der Waal's theory; due to the 
lower abundance and velocities of He atoms this contribution 
is, however, very small and does not influence the calculated 
profiles. Since the theory has not yet been fully extended to 
transitions from ionized species we have to rely on the normal 
van der Waal's broadening approximation by Unsold (1955) 
with an additional enhancement factor E = 2.0 for the Fe II 
lines, which is typical to those adopted in earlier investigations 
of solar Fe II lines (Holweger et al. 1990; Biemont et al. 1991; 
Hannaford et al. 1992). Fortunately, since the Fe II lines are all 
weak, the impact of different choices of E is minor: adopting 
E=1.5 instead leads to only a 0.01 dex increase in the mean 
abundance. Radiative damping was included either with values 
obtained from VALD (Kupka et al. 1999) or calculated from 
the classical formula using the g /-value of the transition; only 
in a few cases does the exact choice of the radiative damping 
influence the results by more than 0.01 dex. Stark broadening 
was not considered. 

A summary of the adopted line data for the Fe I and Fe II 
lines is found in Tables [j], || and |[ The central wavelengths 
for the Fe I and Fe II lines were taken from Nave et al. (1994) 
and Johansson (1998, private communication). When deriving 
elemental abundances this choice is of course of minor im- 
portance, though it is crucial to have accurate estimates when 
studying line asymmetries (Paper I). 

One of the novel features of the current analysis is that the 
Fe abundances are derived from a fit of the line profiles rather 
than from equivalent widths as customary done. This is facil- 
itated by the excellent agreement between observed and pre- 
dicted line shapes, including the departures from perfect sym- 
metry when including the effects of Doppler shifts due to the 
convective flows (Paper I). No microturbulent or macroturbu- 



lent velocities therefore enter the spectral synthesis, since the 
self-consistent velocity field of the simulation is taken into ac- 
count properly. Thereby we have managed to remove three of 
the four hotly debated parameters (W\, £turb and E), which 
have been blamed for the discordance in Fe abundance between 
the Oxford and Hannover- Kiel results. For illustrative purposes 
(e.g. Fig. ^ and||) and for the selection of lines only, we have 
used the quoted equivalent widths from the appropriate sources 
in the literature. We emphasize that they are not used when de- 
termining the Fe abundances. 

For the comparison of line shapes the solar FTS disk-center 
intensity atlas by Brault & Neckel (1987) and Neckel (1999) 
has been used due to its superior quality over the older Liege 
atlas by Delbouille et al. (1973) in terms of wavelength calibra- 
tion (Allende Prieto & Garcia Lopez 1998a,b). In a few cases 
the continuum level was renormalized to better trace the local 
continuum around the lines. 



4. Abundance from weak and intermediate strong Fe I 
lines 

The derived Fe abundances obtained from profile fitting of the 
observed weak and intermediate strong Fe I lines are listed in 
Table [l] We emphasize that the abundances have been derived 
without invoking any equivalent widths, microturbulence or 
macroturbulence, leaving the elemental abundance as the only 
free parameter which is determined by the line strength. The 
good agreement between predicted and observed profiles is il- 
lustrated in Fig. [j]; additional examples can be found in Paper I. 
It is interesting to contrast the remarkable consonance achieved 
with 3D hydrodynamical model atmospheres with the results 
from classical ID model atmospheres, with or without macro- 
turbulence (cf. Fig. |2|; Paper I; Anstee et al. 1997); the improve- 
ment is equally obvious and telling. 

As shown in Fig. |J, the individual abundances show no sig- 
nificant dependence on wavelength or excitation potential. We 
note that the claimed anomalous lines with excitation potential 
of 2.2 eV (Blackwell et al. 1995a) no longer exist in our cal- 
culations. Furthermore, there is no need to fine-tune the tem- 
perature structure to remove existing trends with excitation po- 
tential as necessary with the Holweger-Miiller (1974) model 
with which there is about 0. 15 dex difference in abundances be- 
tween low- and high-excitation Fe I lines (Grevesse & Sauval 
1999). There is, however, a slight trend with line strength in 
the sense that the strongest lines (W\ ~ 9 pm) imply about 
0.06 dex higher abundances than the weakest lines. The reason 
for this behaviour will be discussed further in Sect. [7j although 
it has a minor impact on the final abundance estimates, as illus- 
trated below. It is noteworthy that the trend is more pronounced 
for the Blackwell et al. (1995a) sample than for the lines of 
Holweger et al. (1995), which may partly explain why different 
microturbulences were adopted in the two studies. In terms of a 
ID analysis, the trend would correspond to an underestimated 
microturbulence of about 0.15 km s -1 , which emphasizes the 
relatively minor magnitude of this shortcoming. 
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Fig. 1. A few comparisons between the predicted (diamonds) and observed (solid lines) spatially and temporally averaged Fe I 
lines at disk-center (/i = 1.0). Only every other point in the theoretical profiles are shown for clarity. To illustrate the vast 
improvement over classical ID model atmospheres, in the case of the Fel 615.2 nm line (lower right panel) the corresponding 
prediction with the Holweger-Miiller (1974) model atmosphere (dashed line) when adopting a microturbulence of 0.845 kms -1 
and a Gaussian macroturbulence of 2.4km (the radial-tangential macroturbulence broadening is of course not applicable for 
intensity spectra) is also shown. The Fe abundance for the ID profile has been adjusted to return the same equivalent width as the 
3D profile (log ep c i = 7.59) and the macroturbulence was determined by having the same line depths in ID and 3D. The lack of 
line shift and asymmetry for the theoretical ID profile is clearly seen. Note that all profiles are shown on an absolute wavelength 
scale with no arbitrary wavelength shifts 



The resulting (unweighted) mean abundances of the Oxford 
and Hannover-Kiel samples of weak and intermediate strong 
Fe I lines are log £f c i = 7.46 ±0.04 and log e Fo i = 7.43 ± 0.05, 
respectively, where the quoted uncertainty is the standard de- 
viation (twice the standard deviation of the mean = 0.02). 
The difference between the two samples essentially reflects the 
0.03 dex offset in the absolute scales. From the combined sam- 
ple the estimate is logeFci = 7.44 ± 0.05; it should be noted 
here that nine lines in common have entered twice into this 
result. For reasons outlined above, we consider the transition 
probabilities of Blackwell et al. (1995a) to be slightly superior. 
However, due to the slight trend with line strength our final 
(unweighted) Fe determination is still the mean of all lines: 

loge Fc i = 7.44 ±0.04. 

It should be noted that due to the inclusion of two different 
scales for the oscillator strengths, the scatter is slightly in- 



creased. The final uncertainty is most likely dominated by sys- 
tematic rather than statistical errors, in particular the transition 
probabilities. Furthermore, the neglect of NLTE effects and ob- 
servational complications such as blends and continuum level 
placement, may introduce additional abundance errors which 
are of comparable size; unfortunately astronomy has not yet 
reached the era with 0.02 dex accuracy in absolute abundances, 
in particular not with classical ID model atmospheres, even if 
it is occasionally claimed in the literature. 

The presence of a trend in derived abundances with line 
strength have a minor influence on the mean Fe abundance. 
Restricting the analysis to lines with W\ < 5 pm decreases the 
mean abundances with 0.03 and 0.01 dex for the Oxford and 
Hannover-Kiel samples, respectively, while leaving the scat- 
ter practically intact. The larger sensitivity for the former lines 
can be traced to their in general larger line strengths. This is 
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likely also the reason for the slightly larger difference in mean 
abundances between the two compilations than accounted for 
by the two g/-scales. It should be noted, however, that the Ox- 
ford sample only contains seven lines with W\ < 5 pm, which 
may skew the results somewhat; additional weak Fe I lines with 
high-precision furnace oscillator strengths similar to the pub- 
lished Oxford data would certainly be of great value. 

Given the excellent agreement between the predicted line 
shapes and observed profiles illustrated in Fig. ^ and Paper I, 
very similar abundances to those presented in Table [j] would 
be derived if equivalent widths or line depths had been used 
instead of profile fitting. Due to the larger uncertainties intro- 
duced by the subjectivity of equivalent width measurements 
and departures from LTE in the line cores, such abundance di- 
agnostics are significantly more inferior compared to profile 
fitting, provided of course that the model atmosphere is suffi- 
ciently realistic to accurately predict the line profiles (Asplund 
et al. 2000a). It is interesting to note though that adopting the 
published equivalent widths of Blackwell et al. (1995) would 
result in a 0.02 dex higher Fe I abundance while using the Hol- 
weger et al. (1995) values would result in a 0.03 dex lower 
abundance than those derived from profile fitting for the two 
samples of lines. 

It should be borne in mind that the analysis presented here 
assumes LTE, whose validity may be questioned in particu- 
lar for Fe I lines. Unfortunately no detailed 3D NLTE calcu- 
lations exist for solar Fe lines, and it is therefore difficult to 
predict how the abundances in Table [j] would be altered if de- 
partures from LTE would be allowed. Some preliminary guid- 
ance may come from ID NLTE calculations (e.g. Solanki & 
Steenbock 1988). The calculations by Shchukina (2000, pri- 
vate communication) predict an over-ionization of Fe I and thus 
that the derived ID LTE abundances are slightly underesti- 
mated by 5s 0.1 dex with the Holweger-Miiller (1974) model. 
Simply adopting these ID corrections to our 3D LTE results 
would result in log eFci = 7.50 ± 0.07 for the combined sam- 
ple of Oxford and Kiel lines. Furthermore, the trend with line 
strength would become more pronounced (~ 0.13 dex differ- 
ence between the strongest and weakest lines). Additionally a 
minor trend with excitation potential (~ 0.03 dex difference 
between Xcxc = and 4.5 eV transitions with low-excitation 
lines returning higher abundances) would appear. However, we 
are very reluctant to adopt these results here since it is pre- 
mature to extrapolate ID predictions to the 3D case until 3D 
NLTE calculations for Fe exist. Because the departures from 
LTE depend sensitively on the adopted model atmosphere, the 
temperature inhomogeneities may both amplify or attenuate the 
ID NLTE effects. Naturally such 3D calculations would be of 
great interest. 

5. Abundance from strong Fe I lines 

Strong lines have since long been considered less than ideal 
for the purposes of abundance determinations due to the poorly 
understood collisional broadening which normally requires ad- 
ditional enhancement factors over the classical Unsold (1955) 
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Fig. 2. The derived Fe abundance from weak and intermediate- 
strong Fel lines as functions of wavelength {Upper panel), 
excitation energy (Middle panel) of the lower level, and line 
strength {Lower panel). The lines from the Blackwell et al. 
(1995a) and Holweger et al. (1995) samples are marked with 
• and o, respectively. The dashed lines are linear least square 
fits to the data when including all lines 



recipe. Recent progress in the quantum mechanical treatment 
of the broadening (Anstee & O'Mara 1991, 1995; Barklem 
& O'Mara 1997; Barklem et al. 1998) has, however, opened 
up the possibility to use the damping wings of strong lines, 
which are little sensitive to the non-thermal broadening af- 
fecting weaker lines, as a complement to analyses of weaker 
lines when deriving elemental abundances. Anstee et al. (1997) 
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Fig. 3. A few examples of the predicted (diamonds) and ob- 
served (solid lines) strong Fel lines: Fel 489.07+489.15 (Up- 
per panel) and 523.9 nm {Lower panel) lines. Blending lines 
other than Fe lines have not been included in the synthesis 

found an excellent agreement with the meteoritic abundance 
for the thus determined Fe abundance from strong Fe I lines 
and the Holweger-Muller (1974) model. 

Table || lists the derived Fe abundances with our 3D hydro- 
dynamical solar atmosphere model using a sample of strong 
Fel lines which have been considered the most suitable by 
Anstee et al. (1997) (lines denoted by quality category A+, 
A and A- in their Table 1). Examples of the obtained agree- 
ment between predicted and observed profiles are found in Fig. 
|[ The resulting (weighted) mean Fe abundance from the 14 
strong Fe I lines is 

loge FcI = 7.42 ±0.03. 

However, in spite of being considered as very accurate abun- 
dance diagnostics by Anstee et al. (1997), several of the lines 
turned out to unsuitable due to uncertainties introduced by se- 
vere blending, continuum placement, radiative broadening and 
poorly developed damping wings; those lines are marked in Ta- 
ble ^ and given half weight in the final abundance determina- 
tion. 

Even if the scatter is small for the sample of strong lines, it 
is noteworthy that the standard deviation (er = 0.03) is signifi- 



cantly larger than the claimed accuracy (a = 0.01) of the anal- 
ysis by Anstee et al. (1997) using the Holweger-Muller (1974) 
model. In order to better understand the differences, we have 
therefore re-derived Fe abundances for the same lines in an 
identical procedure to that of Anstee et al. (1997), in particu- 
lar using their collisional broadening data which differ slightly 
from those adopted in Table [| which have been provided by 
Barklem (1999, private communication) from line-by-line cal- 
culations. Our results with the Holweger-Muller (1974) model 
have a significantly larger scatter than that quoted by Anstee et 
al. (1997), identical to our full 3D analysis of the same lines. 
We therefore suspect that the claimed uncertainty in Anstee et 
al. (1997) is over-optimistic and that the true scatter is larger, 
which is also verified by independent calculations by Barklem 
(1999, private communication). The choice of solar atlas (we 
adopt the more recent Brault & Neckel FTS-atlas while Anstee 
et al. use the older Liege atlas) has a minor influence on the 
resulting scatter, although strong lines are often conspicously 
asymmetric in the Liege-atlas (e.g. Ha), presumably due to in- 
accurate continuum tracement. Furthermore there is a system- 
atic offset in theoretical line strengths which amounts to about 
0.03 dex in abundance between our calculations and the identi- 
cal ones by Anstee et al. (J. O'Mara, 1999, private communi- 
cation). The reason for this discrepancy is likely due to slight 
differences in adopted continuum opacities, the P c — P gas - 
relation in the Holweger-Muller (1974) model and code im- 
plementation. This emphasizes again that derived abundances 
rarely have systematic errors smaller than 0.02 dex. 

Due to the subjectivity involved with strong lines in terms 
of choice of solar atlas, continuum placement, wavelength 
shifts, blends, exactly which part of the wings are given the 
greatest weight, and remaining uncertainties in the collisional 
broadening, we consider abundances derived from strong lines 
to be inferior to those from weaker lines, although they natu- 
rally serve as important complements. In this respect it is re- 
assuring that the here derived Fe abundance from strong Fe I 
lines agree well with those from weak and intermediate strong 
Fe I and Fe II lines presented in Sects. ^| and |[ 

6. Abundance from Fe n lines 

The Fe abundance values obtained from the individual Fe II 
lines are listed in Table |J a few examples of the achieved 
agreement between theory and observations are given in Fig. [|. 
As shown in Fig. ^|, the individual abundances show no signifi- 
cant dependence on neither the wavelength, excitation potential 
of the lower level (though the adopted lines provide only a re- 
stricted range), nor the line strength. In this respect, the Fe II 
lines differ from the Fe I lines, which show a minor trend with 
the line strength, at least within the assumption of LTE. 
The resulting (unweighted) mean abundance becomes 

loge FcII = 7.45 ±0.10, 

where the quoted uncertainty is the standard deviation (twice 
the standard deviation of the mean = 0.05). The quoted error 
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lines at disk-center (/i = 1.0) 



of course only reflects the internal accuracy and thereby possi- 
ble uncertainties in e.g. the absolute scale of the ^/-values are 
not accounted for. With Holweger et al.'s ^/-values the mean 
abundance would be 0.04 dex higher while it would be 0.02 dex 
lower with the measurements given in Schnabel et al. (1999); 
the estimated error would remain basically unaltered by such 
exercises. As previously noted, with an Unsold enhancement 
factor of 1.5 instead of 2.0 the abundance would be 0.01 dex 
higher. When excluding the discrepant Fe II 744.93 nm line, 
which is significantly blended in the red wing, the mean abun- 
dance is increased by 0.01 dex; we also note that the <?/-value 
for this line has a comparatively large uncertainty (Hannaford 
et al. 1992). Restricting to the ten lines with W\ < 5.0pm, in- 
creases the mean abundance by only 0.006 dex, demonstrating 
that the observed trend with line strength only affects the Fe I 
lines. 

The main advantage with Fe II lines is their low sensi- 
tivity to details of the temperature structures and departures 
from LTE due to over-ionization. Furthermore, their weakness 
ensures that the lines are formed in the deeper layers which 
are less susceptible to NLTE excitation effects such as photon 
pumping and suction. The abundance derived from Fe II lines 
should therefore be an accurate measure of the solar Fe abun- 



dance, provided the transition probabilities are reliable enough. 
It is reassuring that our average is in good agreement with the 
meteoritic value 7.50 ± 0.01 (Grevesse & Sauval 1998), in par- 
ticular in view of the uncertainty in the absolute (//-scales (Hol- 
weger et al. 1990; Hannaford et al. 1992; Schnabel et al. 1999) 
and that the meteoritic abundance scale probably needs to be 
adjusted downward by about 0.04 dex due to the revised photo- 
spheric Si abundance (Paper III). 

7. Discussion 

The results presented in Sects, ^j, || and ||paint a consistent pic- 
ture for the solar photospheric Fe abundance. Both the weak 
Fe I and Fe II lines suggest very similar abundances: log £f c i = 
7.44±0.05andlogeFo = 7.45±0.10. Since this result does not 
rely on equivalent widths, microturbulence, macroturbulence, 
or, at least for the Fe I lines, collisional damping enhancement 
factors, and is based on highly realistic 3D, hydrodynamical 
model atmospheres, it seems like the long-standing solar Fe 
problem (e.g. Blackwell et al. 1995a,b; Holweger et al. 1995) 
has finally been settled in favour of the meteoritic value, in 
particular considering the slight revision recently of the pho- 
tospheric Si abundance and thus the whole absolute scale for 
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the meteoritic abundances (Paper III): logepo = 7.46 ± 0.01. 
In fact the agreement between the photospheric and meteoritic 
values is partly fortuitious since the remaining uncertainties in 
oscillator strengths and model atmospheres are likely on the 
order of 0.04 dex. 

The difference of 0.1 8 dex (7.64 vs 7.46) between our re- 
sult and the one by Blackwell et al. (1995a) using the same 
set of ^/-values is attributable to a switch to line profile fitting 
and improved collisional broadening treatment, and exchange 
of the microturbulence concept for self-consistent Doppler 
broadening from convective motions and the Holweger-Miiller 
(1974) model for an ab initio 3D hydrodynamical model atmo- 
sphere. Our (unweighted) Fell result logepoii = 7.45 ± 0.05 
is similar to the (weighted) mean logeFoii = 7.47 ± 0.04, 
found by Hannaford et al. (1991) using the same ^/-values 
and the Holweger-Miiller (1974) semi-empirical model atmo- 
sphere, which reflects the small sensitivity of the Fe II lines to 
the details of the model atmospheres; indeed when instead of 
profile fitting the equivalent widths of Hannaford et al. (1991) 
are adopted the (unweighted) mean with the 3D solar model at- 
mosphere is log eFoii = 7.47± 0.04. The main systematic error 
is no longer the model atmospheres and analysis as such, but is 
likely dominated by the accuracy of the transition probabili- 
ties, which still is on the level of 0.03 dex on average for both 
Fe I and Fe II lines, even though the internal precision may be 
higher. 

It is of interest to compare our findings for the solar Fe 
abundance with previously published studies based on 2D 
and 3D hydrodynamical models of the solar photosphere. 
Atroschenko & Gadun (1994) discuss derived Fe abundances 
from Fe I and Fe II lines based on two different types of 
3D model atmospheres (with 30 3 and 32 3 gridpoints, respec- 
tively, to compare with our simulation with the dimension 
200 x 200 x 82) but obtain significantly more discrepant results 
than those presented here: logepd = 7.05 ± 0.06, logeFoii = 
7.48±0.03andloge Fc i = 7.61±0.02, loge Fc „ = 7.42±0.02, 
respectively; here the astrophysically determined ^/-values for 
the Fell lines (using logepo = 7.64) have been rescaled to 
agree with the ones by Hannaford et al. (1992) which we 
have adopted. These results are, however, based on equivalent 
widths for selected lines and the use of microturbulence, which 
had to be introduced in an attempt to hide a very conspicious 
trend with equivalent width. Furthermore, the estimated abun- 
dances only made use of the very weakest lines and therefore 
represent underestimates for the Fe I lines. The discrepancies 
can likely be attributed to the use of gray opacities for the 3D 
model atmospheres and too small height extension, resolution 
and temporal sampling (the spectral synthesis was restricted to 
only one respectively two snapshots from the two simulation 
sequences and therefore should not be considered as proper 
temporal averages). These problems with not sufficiently re- 
alistic model atmospheres are also manifested in the relatively 
poor agreement with observed line profiles and asymmetries. 

The study of Gadun & Pavlenko (1997) suffer from similar 
problems. Their model atmospheres were 2D solar convection 
simulations (with 112x58 gridpoints) but properly temporally 
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Fig. 5. The derived Fe abundance from weak and intermediate- 
strong Fell lines as functions of wavelength {Upper panel), 
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averaged. Utilizing equivalent widths, they derive logeFoi = 
7.33 ± 0.06 and log e Fo ii = 7.44 ± 0.02 for their most reliable 
simulation sequence; again the Fe II result have been rescaled 
for consistency with our analysis. Unfortunately, they do not 
show any comparison between predicted and observed line pro- 
files, but judging from the differences in Fe abundances when 
derived from equivalent widths and line depths we conclude 
that the theoretical line profiles are too narrow, a common prob- 
lem with a too poor numerical resolution in the simulations 
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(Asplund et al. 2000a). To summarize, we are confident that 
our analysis is superior to previously published studies with 
multi-dimensional model atmospheres, a conclusion which is 
further supported by the excellent agreement between the pre- 
dicted line profiles and asymmetries with observations, as de- 
scribed in detail in Paper I, and the confluence between the Fe I, 
Fe II and meteoritic results (Paper III). 

As noted in Sects. ^ and |[ neither the Fe I nor the Fe II re- 
sults depend on the wavelengths or the excitation potentials of 
the lines. Furthermore, the Fe II lines show no trend with line 
strength, in spite of no microturbulence has entered the analy- 
sis, which, as explained in Paper I, is a consequence of the non- 
thermal Doppler broadening from the self-consistently calcu- 
lated convective velocity field. However, according to Fig. |^ the 
individual Fe I abundances appear to depend slightly on the line 
strength, which could signal an underestimated rms vertical ve- 
locity in the line forming layers of the solar simulation (Paper 
I) or a too poor numerical resolution (Asplund et al. 2000a). 
But considering the good overall agreement for the theoretical 
and observed line shapes, which if anything suggests a slightly 
overestimated rms velocity (Paper I) and that no corresponding 
trend is present for the Fe II (Fig. ||) and Si I (Paper III) lines, 
this conclusion seems less likely. Instead we suggest the exis- 
tence of minor departures from LTE in the stronger lines, which 
causes the Fe abundances of these to be slightly overestimated. 
Such departures are more likely to affect Fe I than Fe II lines 
and furthermore stronger lines are more susceptible than weak 
lines due to the decoupling of the non-local radiation field and 
local kinetic gas temperature in the higher atmospheric layers. 
Clearly an investigation of possible NLTE effects for Fe with 
3D inhomogeneous model atmospheres would be interesting, 
similarly to the recent 3D calculations for Li (Kiselman 1997, 
1998; Asplund & Carlsson 2000). 

8. Conclusions 

The application of ab initio 3D hydrodynamical model atmo- 
spheres of the solar photosphere to the line formation of Fe I 
and Fe II lines has allowed an accurate determination of the so- 
lar photospheric Fe abundance. Since such a procedure does 
not invoke any free adjustable parameters besides the treat- 
ment of the numerical viscosity in the construction of the 3D, 
time-dependent, inhomogeneous model atmosphere and the el- 
emental abundance in the 3D spectral synthesis, and consid- 
ering that whole line profiles are fitted rather than equivalent 
widths, the results should provide a more secure abundance 
determination than previously accomplished. The confusion 
introduced by the various choices of mixing length parame- 
ters, microturbulence and macroturbulence no longer needs to 
cloud the conclusions. Furthermore, the analysis has made use 
of recent quantum mechanical calculations for the collisional 
broadening of the Fel lines (Anstee & O'Mara 1991, 1995; 
Barklem & O'Mara 1997; Barklem et al. 1998), which removes 
the problematical damping enhancement parameters normally 
employed, at least for the Fe I lines. In view of these improve- 
ments, it is a significant accomplishment that a consistent pic- 



ture is emerging in terms of Fe abundances: Fe I and Fe II lines 
suggest loge Fe = 7.44 ± 0.05 and loge Fc = 7.45 ± 0.10, 
respectively, which agree very well with the meteoritic value 
logepd = 7.46 ± 0.01 (Paper III) given the remaining un- 
certainties in the transition probabilities. Fe sc i lines may be 
slightly more susceptible for departures from LTE but on the 
other hand the g /-values for Fe II lines are somewhat less ac- 
curate. Our final best estimate for the photospheric Fe abun- 
dance is therefore simply the average of the two results, until 
detailed 3D NLTE calculations and improved measurements of 
the transition probabilities have been performed. Finally, the 
debate of the photospheric Fe abundance (e.g. Blackwell et al. 
1995a,b; Holweger et al. 1995) seems to have been settled in 
favour of the low meteoritic abundance. Also strong Fe I lines 
imply a similar photospheric abundance: logepd = 7.42±0.03 
although we give this result lower weight due to the difficulties 
involved in analysing the wings of strong lines. 

When comparing our results with other recent investiga- 
tions of the solar Fe abundance, it is natural to ask why our 
study should be preferred. After all, traditional analyses using 
classical ID model atmospheres, such as the Holweger-Miiller 
(1974) model, has long been considered sufficient. However, as 
stricter demands are placed on the results in terms of accuracy, 
an improved analysis is required. Why should one embrace 
the results based on hydrostatic ID model atmospheres, equiv- 
alent widths and ad-hoc broadening through microturbulence 
and macroturbulence, when such ID models are inferior to the 
here presented ab initio 3D hydrodynamical models in terms 
of the observational diagnostics available for the Sun, such as 
granulation topology, velocities and statistics, time-scales and 
length-scales of the convection, continuum intensity brightness 
contrast, detailed spectral line profiles, asymmetries and shifts, 
flux distribution, limb-darkening and H-line profiles (e.g. Stein 
& Nordlund 1998; Asplund et al. 1999b; Paper I)? We leave it 
for the reader to ponder this rhetorical question. 
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Table 1. The adopted line data and individually derived abundances for the weak and intermediate strong Fe I lines 



Wavelength a 


xi a 


loj 


19 f 


ref. 


log7 rad a 


lower 


upper 


W b x 


1°£ 


;£fo 


W 


[eV] 






gf b 




level a 


level 3 


[pm] 






438.92451 


0.052 


-4. 


.583 


B 


4.529 


s 


P 


7.17 


7 


.43 


444.54717 


0.087 


-5. 


.441 


B 


4.529 


s 


P 


3.88 


7 


.42 


524.70503 


0.087 


-4. 


.946 


B 


3.894 


s 


P 


6.58 


7. 


.42 


525.02090 


0.121 


-4. 


.938 


B 


3.643 


s 


P 


6.49 


7. 


.45 


570.15444 


2.559 


-2. 


216 


B 


8.167 


s 


P 


8.51 


7 


.53 


595.66943 


0.859 


-4. 


.605 


B 


4.433 


s 


P 


5.08 


7 


.43 


608.27104 


2.223 


-3. 


.573 


B 


6.886 


s 


P 


3.40 


7 


.42 


613.69946 


2.198 


-2. 


950 


B 


8.217 


s 


P 


6.38 


7. 


.46 


615.16182 


2.176 


-3. 


.299 


B 


8.190 


s 


P 


4.82 


7 


.42 


617.33354 


2.223 


-2. 


.880 


B 


8.223 


s 


P 


6.74 


7. 


.44 


620.03130 


2.608 


-2. 


.437 


B 


8.013 


s 


P 


7.56 


7 


.49 


621.92808 


2.198 


-2. 


.433 


B 


8.190 


s 


P 


9.15 


7 


.45 


626.51338 


2.176 


-2. 


.550 


B 


8.220 


s 


p 


8.68 


7 


.45 


628.06182 


0.859 


-4. 


.387 


B 


4.622 


s 


P 


6.24 


7. 


.46 


629.77930 


2.223 


-2. 


.740 


B 


8.190 


s 


P 


7.53 


7. 


.44 


632.26855 


2.588 


-2. 


.426 


B 


8.009 


s 


P 


7.92 


7. 


.51 


648.18701 


2.279 


-2. 


.984 


B 


8.190 


s 


P 


6.42 


7 


.47 


649.89390 


0.958 


-4. 


.699 


B 


4.638 


s 


P 


4.43 


7 


.43 


657.42285 


0.990 


-5. 


004 


B 


4.529 


s 


P 


2.65 


7 


.38 


659.38706 


2.433 


-2. 


.422 


B 


7.936 


s 


P 


8.64 


7. 


.53 


660.91104 


2.559 


-2. 


.692 


B 


7.905 


s 


P 


6.55 


7. 


.49 


662.50220 


1.011 


-5. 


.336 


B 


4.403 


s 


p 


1.36 


7. 


.36 


675.01523 


2.424 


-2. 


.621 


B 


6.886 


s 


P 


7.58 


7 


.48 


694.52051 


2.424 


-2. 


.482 


B 


7.196 


s 


P 


8.38 


7 


.48 


697.88516 


2.484 


-2. 


.500 


B 


6.886 


s 


P 


8.01 


7 


.49 


772.32080 


2.279 


-3. 


617 


B 


6.848 


s 


P 


3.85 


7. 


.55 


504.42114 


2.851 


-2. 


059 


H 


8.009 


p 


s 


7.50 


7. 


.45 


525.34619 


3.283 


-1 


.573 


H 


7.875 


p 


s 


8.10 


7. 


.40 


532.99893 


4.076 


-1. 


.189 


H 


7.659 


d 


P 


5.60 


7 


.47 


541.27856 


4.434 


-1. 


716 


H 


8.226 


p 


d 


1.78 


7 


.45 


549.18315 


4.186 


-2. 


.188 


H 


8.158 


d 


p 


1.06 


7 


.41 


552.55444 


4.230 


-1. 


.084 


H 


8.382 


P 


s 


5.80 


7. 


.39 


566.13457 


4.284 


-1. 


.756 


H 


7.908 


P 


s 


1.98 


7 


.38 


570.15444 


2.559 


-2. 


.130 


H 


8.167 


s 


P 


8.60 


7. 


.45 


570.54648 


4.301 


-1. 


.355 


H 


8.290 


p 


s 


3.90 


7 


.38 


577.84531 


2.588 


-3. 


440 


H 


8.167 


s 


p 


1.95 


7 


.36 


578.46582 


3.396 


-2. 


.530 


H 


7.877 


p 


s 


2.50 


7 


.39 


585.50767 


4.607 


-1. 


.478 


H 


8.281 


p 


d 


2.10 


7. 


.41 


608.27104 


2.223 


-3. 


590 


H 


6.886 


s 


P 


2.82 


7. 


.44 


615.16182 


2.176 


-3. 


.270 


H 


8.190 


s 


p 


4.56 


7 


.39 


621.92808 


2.198 


-2. 


.422 


H 


8.190 


s 


P 


8.70 


7 


.44 


624.06460 


2.223 


-3. 


.230 


H 


7.196 


s 


P 


4.38 


7 


.36 


627.12788 


3.332 


-2. 


.703 


H 


8.074 


p 


s 


2.09 


7 


.40 


629.77930 


2.223 


-2. 


.727 


H 


8.190 


s 


P 


7.30 


7. 


.42 


648.18701 


2.279 


-2. 


960 


H 


8.190 


s 


P 


6.30 


7 


.45 


658.12100 


1.485 


-4. 


680 


H 


7.193 


s 


p 


1.41 


7. 


.39 


666.77114 


4.584 


-2. 


.112 


H 


8.158 


s 


P 


0.89 


7 


.55 


669.91416 


4.593 


-2. 


.101 


H 


8.158 


s 


p 


0.73 


7 


.45 


673.95220 


1.557 


-4. 


.790 


H 


7.176 


s 


P 


1.03 


7 


.30 


675.01523 


2.424 


-2. 


610 


H 


6.886 


s 


P 


7.70 


7. 


.47 


679.32593 


4.076 


-2. 


.326 


H 


7.622 


d 


P 


1.10 


7. 


.39 


680.42715 


4.584 


-1 


.813 


H 


7.719 


s 


P 


1.40 


7 


.46 


683.70059 


4.593 


-1. 


687 


H 


7.719 


s 


P 


1.54 


7 


.44 


685.48228 


4.593 


-1. 


.926 


H 


7.659 


s 


P 


1.00 


7 


.52 


694.52051 


2.424 


-2. 


440 


H 


7.196 


s 


P 


8.20 


7 


.44 


697.19330 


3.018 


-3. 


.340 


H 


8.161 


s 


P 


1.20 


7. 


.36 


697.88516 


2.484 


-2. 


.480 


H 


6.886 


s 


P 


7.90 


7. 


.47 


718.91510 


3.071 


-2. 


.771 


H 


8.161 


s 


P 


3.80 


7. 


.53 


740.16851 


4.186 


-1. 


599 


H 


7.847 


d 


P 


4.10 


7 


.50 
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Table 2. The adopted line data for the strong Fe I lines 


Wavelength a 
[nm] 


xi a 

[eV] 


log 5/ 


ref. 
gf b 


log7 rad a 


lower 
level 


upper 
level 


cr c 




a c 


log e Fo 


407.17380 


1. 


.608 


-0.022 


B 


8.009 


s 


P 


328 


0. 


.252 


7.40 d 


438.35449 


1 


485 


0.200 


B 


7.936 


s 


P 


295 


0. 


.265 


7.43 d 


441.51226 


1 


.608 


-0.615 


B 


7.986 


s 


P 


305 


0. 


.261 


7.46 d 


489.07549 


2. 


.875 


-0.390 





8.004 


p 


s 


747 


0. 


.236 


7.45 


489.14922 


2. 


.851 


-0.110 


O 


8.009 


p 


s 


739 


0. 


.236 


7.45 


491.89941 


2. 


.865 


-0.340 


O 


8.009 


p 


s 


739 


0. 


.237 


7.46 d 


495.72988 


2. 


.851 


-0.410 


O 


8.009 


p 


s 


727 


0. 


.238 


7.43 


495.75967 


2. 


.808 


0.230 


O 


8.009 


p 


s 


713 


0. 


.238 


7.43 


523.29404 


2 


.940 


-0.060 


O 


8.009 


p 


s 


712 


0. 


.238 


7.41 


526.95376 





.859 


-1.321 


B 


7.185 


s 


p 


237 


0. 


.249 


7.38 


532.80386 





.915 


-1.466 


B 


7.161 


s 


P 


239 


0. 


.248 


7.40 


532.85317 


1 


.557 


-1.850 


O 


6.848 


s 


P 


282 


0. 


.252 


7.40 d 


537.14897 





.958 


-1.645 


B 


7.152 


s 


p 


240 


0. 


.248 


7.41 d 


544.69170 





.990 


-1.910 


O 


7.152 


s 


p 


241 


0. 


.248 


7.40 d 



a From Nave et al. (1994) and the VALD data base (Kupka et al. 1999) 

b From Blackwell et al. (1995a) (ref. gf=B) and O'Brian et al. (1991) (ref. gf=0) 

c Collisional broadening data from Barklem (1999, private communication) 

d Lines which are given half weight in the final abundance estimate due to uncertainties introduced by blending lines, continuum placement, 
radiation broadening and poorly developed damping wings 



